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Abstract: We consider 3+1-dimensional gauge theories at finite temperature and a finite 
density of charges which couple to a 2+1-dimensional Chern-Simons operator, giving rise 
to a #-term with constant spatial gradient of 6. The strong-coupling limit of thermal 
J\f = 4 super- Yang- Mills theory with this kind of anisotropic deformation has been used in 
the context of the AdS/CFT correspondence as a model for strongly coupled anisotropic 
quark-gluon plasma. In this paper we work out the thermodynamics and the (nontrivial) 
phase diagram in the limit of vanishing gauge coupling and compare with the corresponding 
strong-coupling results. 
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1. Introduction 

There is by now little doubt that quark-gluon plasma can be produced and studied in 
current ultrarelativistic heavy-ion collider experiments (RHIC at BNL, LHC at CERN) [1]. 
The quark-gluon plasma created in this way is however initially very far from equilibrium, 
and it may have substantial pressure anisotropics over its entire lifetime until freeze-out 
[2, 3, 4, 5]. 

At weak coupling, an anisotropic plasma exhibits plasma instabilities [6, 7, 8, 9, 10] 
which in the non-Abelian case have complicated dynamics that has been studied extensively 
by numerical approaches using the approximation of stationary anisotropy [11, 12, 13, 14, 
15, 16, 17], and more recently also for anisotropic expansion [18, 19, 20]. 1 These plasma 
instabilities are crucial for understanding thermalization and isotropization of a weakly 

1 Similar instabilities have been identified in the so-called glasma phase, where the dynamics of non- 
Abelian gauge fields after the collision of ultrarelativistic color sources is governed in leading order of 
perturbative QCD by classical Yang-Mills field equations [21, 22] as well as in classical-statistical simulations 
of Yang- Mills field dynamics [23, 24]. 
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coupled plasma [9, 25, 26]. Moreover, they lead to "anomalously" low effective viscosity and 
could thus mimic an inherently strong-coupling situation [27] . It is therefore of considerable 
interest to develop a more complete understanding of the effects of anisotropies in both 
weakly and strongly coupled plasmas. 

At strong coupling, where the AdS/CFT correspondence [28] has provided new tools 
[29], the effects of (temporarily fixed) anisotropy in a supersymmetric Yang-Mills plasma 
have been modeled by singular geometries involving comparatively benign naked singulari- 
ties [30, 31] and more recently by a completely regular construction involving axion-dilaton 
gravity [32, 33]. Interestingly enough, this latter model has been found to contain phases 
with instabilities reminiscent of the filamentation instabilities at weak coupling. In follow- 
up works, several observables of interest to heavy-ion physics have been studied in this 
model [34, 35, 36, 37, 38] and compared with calculations for a weakly coupled anisotropic 
plasma [39, 40, 41, 42, 43, 44, 45], a perhaps particularly remarkable finding being that of 
a shear viscosity coefficient breaking the Kovtun-Son-Starinets bound [46, 47]. 2 

While other holographic models for anisotropic fluids have been constructed [49, 50, 
51], the model of Ref. [32, 33], which builds upon the string theory dual constructed in 
Ref. [52], is especially attractive because it is a string-theoretic top-down construction. 
The stationary anisotropy is brought about by an anisotropic distribution of D7 branes 
which apart from wrapping the S 5 of the bulk geometry fill only 2 out of the 3 spatial 
dimensions. Like the color D3 branes they do not extend along the holographic direction, 
but are dissolved in the geometry. The corresponding gravitational background involves 
an axion with constant spacelike gradient. 

On the gauge theory side, this corresponds to a deformation of N = 4 super- Yang-Mills 
theory by a position-dependent #-term 



5S= I ^TrFAF= / -^TrFAF (1.1) 
^ J S YM 



where z is the spatial coordinate along which a constant density ud-j of D7 branes is set 
up. Absorbing gym in i 7 by a rescaling this leads to the modified Yang-Mills Lagrangian 

C = -\f%F^ - \e{x)e^F« u F a pa (1.2) 

with 9{x) = az. By partial integration one finds that the charge density a is coupled to a 
homogeneous but anisotropic operator given by a 2+1-dimensional Chern-Simons term. 

The aim of the present paper is to study the thermodynamics and phase diagram of 
this anisotropic Chern-Simons deformed gauge theory at vanishing coupling and to compare 
with the corresponding results obtained by means of gauge-gravity duality in Ref. [33] . As 
we shall see, the phase diagram at zero coupling is even richer than the one at strong 
coupling, with significant differences in particular at high temperatures, T a. 

A similar Chern-Simons deformation of electrodynamics has originally been studied 
by Carroll, Field, and Jackiw in [53] as a model for a Lorentz and CPT violating electro- 
dynamics which preserves rotational invariance by a timelike gradient of the 9 parameter. 



2 See [48] for a related finding in a different gauge-string dual. 
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In this version of the theory there are tachyonic modes, which are however absent in our 
case of interest, the case of a spacelike 9 gradient [54, 55] (which for a while has attracted 
attention as a possible explanation for cosmic anisotropy in the polarization of distant radio 
sources [56] but evidence for the latter was refuted by Carroll and Field [57, 58]). 

As we shall discuss further below, the photons in this modified anisotropic electrody- 
namics have dispersion laws of the form 



c4 = k" + y ^l±yi + -/|. (1-3) 

There are gauge boson modes with a mass gap a as well as ungapped ones, but whenever 
there is a wave vector component parallel to the direction of anisotropy, there is a deviation 
from an ordinary mass shell or the light-cone. While ungapped modes have w_ < |k|, all 
modes have u± > and are therefore stable, which is in fact markedly different from the 
gauge boson propagator in a weakly coupled anisotropic plasma, where a rich spectrum of 
instabilities and unstable modes arises [8, 9]. Nevertheless, like in the strong coupling case 
we shall identify phases with thermodynamic instabilities against inhomogeneous redistri- 
bution of the "Chern-Simons charge" density a. 

2. Setup and notation 

We consider a theory of free photons in a system containing a source j(x) for 2+1 di- 
mensional Chern-Simons operator. The partition function is given by the Euclidean path 
integral (in Euclidean spacetime, we write all the Lorentz indices as subscripts to distin- 
guish it from Minkowski spacetime) 



Z(T,j) = J VA^exp - ' dr J d 3 x (c + ^j(x)e pupa A p F up ( a + 



(2.1) 



C = -F pv F pv + - (d p A p ) + ghosts, (2.2) 



l F F + 1 

A r iiv r ii,i> ' 2 

where the constant spacelike unit vector specifies the 2+1 dimensional subspace, the 
factor 1 /4 is a convention, SI is the cosmological constant (needed for renormalization) and 
we have adopted the Feynman gauge. Note that the charges that couple to e^p^A^F^^a- 
can be localized only in the (-direction, j(x) = ftji/A otherwise gauge invariance is 
broken. We consider a system in which we have a constant charge density along the (- 
direction (which without loss of generality can be taken as the z-direction) , 

j( z) = a = (2.3) 
L ll 

so that (2.1) is anisotropic but translationally invariant. Here, N is the number of charges 
coupling to the Chern-Simons operator and Lm the extent of the system in ("-direction. In 
order to make contact with the notation of Ref. [32, 33], we write for the constant density 
j(z) = a. Denoting the linear extent of the system perpendicular to (-direction as L±, the 
volume is then V = L 2 , L\\ . 
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We use the notation K 2 = k^k^ = k^ + k 2 = k 2 + k\ + k 2 for momenta. We renor- 
malize the theory in dimensional regularization by splitting the spacetime into R 3 x R 1_2e 
("transverse" x "longitudinal"), with the Chern-Simons term taking the form e^AiFj^ 
with {i,j,k} labeling the transverse directions (including Euclidean time). Additionally, 
we label the longitudinal directions with letters from the beginning of the Greek alphabet, 
a, f3, . . . , leaving fi,u, . . . to label directions in the entire spacetime. 

Our notation for momentum integrations is such that 

d 2 k ± fe^A 2 Y f d l -^k\\ d 2 k ± 

(2.4) 



i 



(2tt) 2 \ A-k ) J (27T) 1 - 26 (2tt) 2 ' 
TJ2 I (2-5) 



K k - 
3. Thermodynamic functions 

Before carrying out explicit calculations, we discuss the various thermodynamic functions 
of interest to us (see also Appendix C of Ref . [33] ) . 

Energy is a function of extensive variables, E = E(S, N, L±, Ly). Pressures in various 
directions are most naturally expressed in terms of it, and are given by 

p±= -^(S>w " , = "4(^) w / (3 - 1) 

A more convenient quantity to compute, however, is the Helmholtz free energy, F = 
F(T, N,L ± ,L ]] ) = E - TS, given by 

F = -TlnZ = Vf(T,N/L\\) = Vf(T,a). (3.2) 

(The functional form F(T, N, L±, Ln ) = Vf(T, N/Lu ) is clear from the form of the partition 
function.) From the definitions of pressures in Eq. (3.1) and the relation between F and 
E, we immediately find that the pressures are given by 

1 / dF \ 1 / dF 

T.N.Lu L \\ \ dL ±/T,a 



Note that even though our physical system is inherently anisotropic as soon as a ^ 
(with, e.g., anisotropic relations between energies and momenta of its particles, see be- 
low), the pressure would be necessarily isotropic in thermal equilibrium if a did not react 
differently to changes in the system size along different directions. 

We can also study the system in the grand canonical ensemble, as a function of the 
chemical potential /i = dF/dN conjugate to the number of the Chern-Simons charges. 
The associated free energy, the grand potential, is given by 



\8N J t,L\,Lu 



f 



da 



(3.5) 
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The grand potential has the general form G(T, fj,, L±,L») = Vg(T, fi/L 2 ± ), which is easy to 
see from above as follows. Defining 

(3.6) 

oa 

it is clear by construction that the grand potential can be written as 

Of 

G = Vg(T,<S>) with g (T,<f>) = f-a^ a . (3.7) 



On the other hand, we have 



^ V [M) TiL = Li {ta) T = L ^ (38 » 

and thus $ = n/Lj_ and G = Vg(T, fj,/L±). In terms of the grand potential, the pressures 
can then be written as 

Instead of working with F and G, which are functions of conjugate variables N and /x, 
respectively, it is more convenient to consider the free energy densities / and g, which are 
functions of the conjugate variables a (number of charges per unit length) and <3? (chemical 
potential per unit transverse area), respectively. 

4. Results at vanishing gauge coupling 
4.1 Free energy 

Since we are considering a free theory, the path integral in Eq. (2.1) for a homogeneous 
source j(z) = a can be carried out using standard methods. Details of all the computations 
are given in Appendix A; we will quote only the essential results here. The resulting 
Helmholtz free energy density f(T,a) = — T/V In Z(T, a) is given by 



± 

where 



f(T,a) = Sl + TY J j k ^± + ln(l-e 



-Poj± 



(4.1) 



a 2 ± a 4 + 4a 2 k?, 

ul = kl + M£(k\\), M 2 (k {l ) = k 2 + 1 ) (42) 

or, equivalently, 

M ± m) = ^k 2 + a 2 /±± a -. (4.3) 

Note that (4.1) is an even function of a so that without loss of generality we can assume 
a > 0. 
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The T = limit of Eq. (4.1) contains a UV-divergence that needs to be renormalized. 
In the MS-scheme, we obtain 



/(0,o) 



where 



ca 5a , a 

12^ + 256^1 + " (A) ' 



dx (x 2 + i(l + v 7 ! + 4x 2 )) 7 + (V + 1(1 - v 7 ! + 4x 2 )^) 



(4.4) 



\ 3/2 

> 



-2x d 



15 



64V1 + x 2 



0.29136... 



(4.5) 



and fi(A) is the renormalized cosmological constant running with the scale A. Note that 
Eq. (4.4) yields the exact renormalization group equation governing the running of the 
cosmological constant in this model. However, since we are primarily interested in the 
thermodynamics as a function of a at some fixed (but arbitrary) scale A, we are free to 
choose the value of the cosmological constant (which by definition is independent of a) at 
that scale. With this in mind, fixing A to give units a is measured in, we set fi(A) = 0, 
yielding for the T = free energy 



f(0,a) 



ca 5a a 
12vr 2 + 256vr 2 n A' 



The coefficient in front of the logarithm gives the trace anomaly of our system, 

5a 4 



e ~ 2p± - P " = - 25^ ■ 



(4.6) 



(4.7) 



Note that the free energy vanishes at a finite value of the density, 3 a = ao, and we can 
thus express the scale A in terms of ao, 



such that 



In 



f(0,a) 



ao 
A 



64c 
15 



5a a 
2567T 2 ag 



(4.8) 
(4.9) 



From now on, we express everything in units of ao- 

Other thermodynamic variables of interest at T = can now be computed, 



0£ 
da 

d<S> 
da 



5a 3 



256vr 2 
5a 4 

256vr 2 
5a 2 

256vr 2 



1 +41n 
1 + 3 In 



7 + 12 In 



a 



(4.10) 
(4.11) 
(4.12) 



3 In case we had not set the cosmological constant to zero, the equivalent statement is that the free 
energy has the same value at some finite density a — ao as at a — 0. 
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Respectively, they vanish at different densities a given by 



a* = a e" 1/4 « 0.7788 a , (4.13) 

p\\ 



a = a e~ 1/3 « 0.7165 a , (4.14) 



a*/ = a e" 7/12 « 0.5580 a . (4.15) 
The finite-T contribution is given by 

f(T, a) - /(0, a) = TJ2 ^ In (l - e"^±) . (4.16) 

Unfortunately, it is not possible to express this sum-integral in a closed form. We can, how- 
ever, derive a simpler integral representation for the finite-T contribution that is straight- 
forward to evaluate numerically (see Appendix A. 2 for details), 

where y = a/(2T). The integral in (4.17) is of the same form as integrals encountered in 
the thermodynamics of bosons of mass a/2 and chemical potential —a/2 (specifically, the 
first term gives exactly the free energy of such bosons), and we can use standard methods 
to retrieve high- and low temperature expansions of it. At high temperatures (T S> a) we 
obtain 

ir 2 T* a 2 T 2 a 3 T 5a' f, a 1 \ . _, „ 



f(T, a) - f {0 , a) _ -_ + _ - _ - — 2 _ + 7£ - - j + 0(a«), (4.18) 
while at low temperatures (T <C a) we get 

/(T, a) - /(0, a) = a 1 ^ + 0{a~^). (4.19) 

These expansions are useful in studying the asymptotic behaviour of the system, but in the 
following analysis we solve the system numerically using the exact results in Eqs. (4.17) 
and (4.9). 

4.2 Phase diagram and instabilities 

The system we consider contains a fixed number N of charges that couple to the Chern- 
Simons term. However, how these charges are distributed along the z-axis is a free param- 
eter determined by minimizing the free energy F. We refer to these various distributions 
as various phases of the system. For simplicity, we only consider distributions which con- 
sist of homogeneous regions of finite extent along the direction of anisotropy with various 
densities, since we know how to compute the free energy density of such homogeneous 
regions. We refer to the phase with just one region as homogeneous phase, and the phase 
with two or more homogeneous regions as inhomogeneous phase. Note that the boundaries 
between regions in the inhomogeneous phase must be perpendicular to the z-axis, since 
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gauge invariance requires a to be just a function of z. Hence, any resulting inhomogeneous 
structure is "lasagne-like" . 

In cases where the homogeneous phase is energetically disfavored, it can still be 
metastable, which means that finite charge fluctuations are required for the transition 
to the inhomogeneous phase. When arbitrarily small charge fluctuations are sufficient for 
lowering the free energy, the system is called thermodynamically unstable. 

4.2.1 Coexistence of homogeneous regions 

In the inhomogeneous phase, the free energy can be minimized by a redistribution of the 
given number of charges into an arbitrary number of homogeneous regions. We can analyse 
this situation by considering the case of two neighboring homogeneous regions. 

Two separately homogeneous regions can coexist if the intensive thermodynamic vari- 
ables have the same values in each region. In the present case this means regions with 
different values of a can coexist if the longitudinal pressure py has the same value in each 
region ("mechanical equilibrium"). Since the boundary between regions must be perpen- 
dicular to the z-axis, we do not need to require equality of the transverse pressures. In 
addition, also the chemical potential $ must be the same in each region ("chemical equi- 
librium"), as long as each region contains a non-zero number of charges so that processes 
taking charges from one region to another can remain in equilibrium. If one region is empty 
of charges (and hence isotropic), chemical equilibrium cannot be reached. 4 




Figure 1: Longitudinal pressure as a function of 1/a for various temperatures, normalized to 
Stefan-Boltzmann pressure. 

In Fig. 1 we have plotted the longitudinal pressure as a function of 1/a (= Ly/iV, 
so this corresponds to the traditional p-V diagram) for various temperatures. As can be 
seen, for sufficiently low temperatures, there are various values of a that yield the same 

4 Charges that couple to the Chern-Simons term are not dynamical in the model, which means that 
formation of charge-anticharge pairs is excluded. Hence, regions with negative charge density cannot emerge. 
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longitudinal pressure. Moreover, it can immediately be seen that the homogeneous phase 
is thermodynamically unstable for certain values of a, since the condition 



< 



( d n 

\d(l/a) 



< o, 



^ dL \\J T ,N 

would be violated. Equivalently, we can express this condition as 



(4.20) 



da 2 J T 



( Of! ) rp > 



(4.21) 



A homogeneous phase violating this condition cannot exist in equilibrium. We can thus 
conclude that, for sufficiently low temperatures and certain values of the overall density a, 
the system must reside in the inhomogeneous phase. On the other hand, for sufficiently 
large temperatures, the system must reside in the homogeneous phase since the longitudinal 
pressure is then a monotonically increasing function of a. 

More precisely, the conditions for coexistence of two homogeneous regions with densi- 
ties a± and a 2 , < a\ < a 2 , are given by: 

(a) coexistence of two anisotropic regions (a>i > 0): 



$(T,ai) = $(T,a 2 ) = $ 

/(r,a 2 )-/(r, Ql ) = g 

a 2 — a\ 

(b) coexistence of an isotropic and an anisotropic region (a\ = 0): 



(4.22) 
(4.23) 



f(T,a 2 )-a2$(T,a 2 ) = f(T,0) 



7T 2 T 

45 



(4.24) 



/ 



/ 



Homogeneous 




Inhomogeneous 










a 2 





Homogeneous 
Inhomogeneous 


«2 / 


"~ ■«. ^ 





(a) 



(b) 



Figure 2: Two cases where the conditions for coexistence of homogeneous regions are realized, 
(a) The dashed line, giving the free energy of the inhomogeneous phase, is the tangent to f(T,a) 
at points a = a\ and a — a<i- (b) The dashed line is the tangent to f(T, a) at a = ai and coincides 
with f(T,a) at a = 0. 
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These two cases are graphically represented in Fig. 2. Two anisotropic regions with 
densities a± and 02 can coexist if the straight line drawn through points (a±, f(T, a%)) and 
(ct2j f(T, 02)) is the tangent to the f(T,a) curve at the points a± and a<i (Fig. 2a). That 
straight line corresponds to the free energy of a system containing a mixture of homogeneous 
phases with densities a\ and 02, that is, the free energy of the inhomogeneous phase. If 
that line lies below the f(T,a) curve, then the free energy of the inhomogeneous phase is 
lower than that of homogeneous phase with corresponding overall density, and thus the 
equilibrium state is given by the inhomogeneous phase. Since stability of the homogeneous 
regions requires that d 2 f /da 2 \ ai ,a 2 > 0, at least part of the f(T,a) curve must lie above 
the line. On the other hand (Fig. 2b), an isotropic region void of the Chern-Simons charges 
can coexist with an anisotropic region with charge density a = 02 if the line drawn from 
(0, /(T, 0)) to (02, /(T, 02)) is tangent to f(T,a) at a = 0,2- This guarantees mechanical 
equilibrium. 

4.2.2 T-a phase diagram 
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Figure 3: The phase diagram. The solid line separates the homogeneous and inhomogeneous 
(shaded region) phases. In the green shaded part of the phase diagram, the inhomogeneous phase 
consists of separate anisotropic regions with different values of a, whereas in the red shaded part, the 
plasma contains anisotropic and isotropic regions. The red dashed line indicates the region where the 
(would-be) homogeneous phase is unstable. The blue dash-dotted line indicates vanishing chemical 
potential. Inside the line, the pressure anisotropy is oblate, outside it is prolate. The dotted line 
indicates the region where f(a,T) < f(0,T). 



The phase diagram of the stable phases can now be computed, and is given in Fig. 3. 
The shaded area indicates the region in phase space where thermodynamic equilibrium is 
reached in the inhomogeneous phase with a mix of values a\ and 02 given by the boundary 
values of a of the shaded region at the given temperature. Within the inhomogeneous 
region, we have indicated with a dashed line the region where the homogeneous phase 
would be thermodynamically unstable. Outside the dashed line the homogeneous phase is 
at least metastable. In the red shaded area, the inhomogeneous phase consists of isotropic 
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and anisotropic (finite Chern-Simons charge density) regions. In the green shaded region 
the system is, in addition to being inhomogeneous, also anisotropic everywhere. 

The type of (pressure) anisotropy of the system (prolate vs. oblate) is determined by 
the sign of the chemical potential & = df /da. For $ > 0, the plasma is prolate, i.e, 
P\\ > P±: see Eq. (3.4). In Fig. 3, the boundary in phase space separating prolate plasma 
from oblate plasma is given by the blue dot-dashed line. Outside that line, the plasma is 
prolate. Note that in the red shaded part of the phase diagram, the chemical potential does 
not have a unique value within the plasma: the plasma is composed of isotropic ($ = 0) 
and oblate (<I> < 0) regions. In the green shaded part of the phase diagram, the plasma 
contains only prolate ($ > 0) regions, which are in chemical equilibrium (same <!> but 
different a). 

In the homogeneous phase (unshaded region in Fig. 3), the plasma is prolate for the 
major part of the phase diagram. There is only a narrow region at low temperatures where 
the plasma is homogeneous and oblate (the unshaded region between the full line and the 
blue dash-dotted line). 

Note that we have chosen to characterise the plasma as oblate, prolate or isotropic only 
with regard to the ratio of longitudinal and transverse pressures. Intrinsically the plasma 
is anisotropic whenever a / 0, but it can have isotropic pressure also for nonvanishing a, 
when <I> = (along the blue dash-dotted line in Fig. 3). 

4.2.3 Metastable homogeneous phases 




Figure 4: The phase diagram showing different metastable phases labelled A-D with homogeneous 
Chern-Simons charge distribution (the thermodynamically unstable region is shown in yellow). 
The light-blue dash-dotted curve corresponds to $ = in the homogeneous case, with $ > 
(prolate pressure anisotropy) and $ < (oblate) above and below this line, respectively. Above 
the straight dark-blue dash-dotted line the inhomogeneous ground state is composed of differently 
prolate plasma, below this line the ground state is a mix of isotropic and oblate plasma. 



In the region of the phase diagram where the inhomogeneous phase is energetically 
preferred (in the area below the solid line in Fig. 3), one can distinguish a number of 
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qualitatively different metastable situations for a homogeneous charge distribution as shown 
in Fig. 4. In this diagram, the metastable homogeneous phases are labelled A-D while the 
region where the homogeneous phase is thermodynamically unstable is colored yellow. For 
a homogeneous charge distribution the dividing line between oblate and prolate pressure 
anisotropy is given by light-blue dash-dotted curve, whereas for the inhomogenous phase 
the dividing line is the straight dark-blue dash-dotted line. This difference is responsible 
for the appearance of four different "decay modes" of homogeneous metastable phases with 
density a: 

(A) the metastable homogeneous phase is oblate and decays into a mix of regions that 
are isotropic (ai = 0) and oblate (02 > a) 

(B) the metastable homogeneous phase is oblate and decays into a mix of prolate regions 
with different nonvanishing densities a\ < a and 02 > a 

(C) the metastable homogeneous phase is prolate and decays into a mix of prolate regions 
with different nonvanishing densities a\ < a and 02 > a 

(D) the metastable homogeneous phase is prolate and decays into a mix of regions that 
are isotropic (ai = 0) and oblate (02 > a) 

Note that the metastable homogeneous phases of type A and D, which like all homoge- 
neous phases are trivially in chemical equilibrium, decay into inhomogeneous systems that 
are no longer in chemical equilibrium. 

5. Comparison with holographic infinite-coupling results 

In Ref. [32, 33] the gravity dual to the strong-coupling limit of maximally supersymmetric 
Yang-Mills theory at infinite color number N c with Chern-Simons charge density a has 
been constructed and worked out in great detail. 

In the 5-dimensional bulk, this involves a linear axion field x = az an d an anisotropic 
metric of the form 

ds 2 = - F(u)B(u)dt 2 + dx 2 + dy 2 + U(u)dz 2 + — - , (5.1) 

which at finite temperature has a regular horizon at some value u = u^. and approaches 
anti-de Sitter form at the boundary u = 0. 

At nonvanishing a, the gauge theory has a trace anomaly, which divided by N 2 
reads [33] 

(e-2p ± - Pll )/N 2 = -^- 2 . (5.2) 

This is in fact curiously close to the zero-coupling (and nonsupersymmetric) result (4.7) 
which equals — a 4 /(51.2 tt 2 ). 

When expressed in terms of ao, defined such that /(0, ao) = at nonvanishing a's, 
the result for / at strong coupling is just proportional to the zero-coupling case, which 
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implies that a$, a p ,, , a$' are exactly the same as in Eqs. (4.13)-(4.15). Therefore, the phase 
diagrams of the two theories when drawn in units of ao exactly coincide in the limit of 
vanishing temperature. 

At nonzero temperature, there are however significant differences. At low temperatures 
(T <C a), the free energy behaves as [33] 

[f(T,a) - f(0,a)]/N 2 ~ -0.9 a^T 11 / 3 (5.3) 

while at zero coupling we had obtained -SmQa 1 / 2 ^/ 2 in Eq. (4.19). 
At high temperatures, the strong coupling result reads [33] 

[/(T,a) - f(0,a)]/N 2 = - ^ + 0(a 4 ), (5.4) 

which is to be contrasted with Eq. (4.18), where we had -7r 2 T 4 /45 + a 2 T 2 /48 - a 3 T/64 + 
0(a 4 ). 

The term cubic in a mass parameter and linear in T in the zero-coupling result is 
typical of a weakly coupled plasma involving massive quasiparticles; its absence in the 
strong-coupling result illustrates that a quasiparticle description is no longer possible there. 

The different sign in the a 2 T 2 correction implies a different sign for $ = df/da, which 
is solely responsible for the pressure anisotropy. This means that in the high-temperature 
limit the strong-coupling system has oblate pressure anisotropy, whereas it was prolate at 
zero coupling. This also turns out to completely change the structure of the phase diagram 
at high temperature. 

In Fig. 5 we have evaluated numerically the longitudinal pressure of the strong-coupling 
theory as a function of 1/a for the same set of temperatures (in terms of the scale ao) as 
in the zero-coupling theory (Fig. 1). In contrast to the latter we now find that increasing 
the temperature does not get rid of unstable regions in a; the homogeneous phase is now 
always unstable for sufficiently small a while stable for larger a. 
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Figure 7: The phase diagram at strong coupling. Labeling of the different curves as in Fig. 3. 



Fig. 6 shows a typical dependence of f(a, T) on a. Because the curvature of / is 
negative at a = 0, but / behaves as a 4 In a for sufficiently large a, there is always a region 
in a around a = where the free energy can be lowered by an inhomogeneous mix of 
isotropic (a = 0) and more strongly anisotropic (a = 02) domains. 

In Fig. 7 our numerical result for the full phase diagram at strong coupling is shown, 
which agrees with the qualitative sketch given in Ref. [33] 5 . Comparing with the zero- 
coupling result shown in Fig. 3 we see that at small temperatures the two phase diagrams 
are very similar, but at high temperatures there are significant differences. 

In the strong-coupling case, instabilities of homogeneous phases against the formation 
of inhomogeneous structures along the z-direction only occur in the region left to the 
blue dash-dotted line where the pressure anisotropy is oblate, but the oblate region now 
extends to arbitrarily high temperature. Between the blue dash-dotted line and the full 

5 The particular bending of the curves at low temperatures has been greatly exaggerated in Ref. [33], 
but is qualitatively correct apart from the fact that all lines become strictly vertical in the limit T — > 0. 



- 14 - 



line the oblate homogeneous plasma is stable. In the red-shaded region the energetically 
preferred inhomogeneous phase always corresponds to a mixture of isotropic (a = 0) and 
more strongly anisotropic regions. To the left of the red dahed line, the homogeneous phase 
is thermodynamically unstable, to the right it is at least metastable. Only type A of the 
metastable phases discussed in Sect. 4.2.3 is realized at strong coupling. 

As we have seen above, in the weak-coupling case, sufficiently high temperature implies 
prolate pressure anisotropy and stability of a homogeneous distribution of Chern-Simons 
charge. However, at moderate temperature, there are also domains with prolate anisotropy 
that are unstable as well as inhomogeneous phases with a mixture of prolate plasma with 
different nonzero Chern-Simons charge densities (the green-shaded area in Fig. 3), which 
do not occur in the strong-coupling case. 

6. Conclusion 

We have worked out the complete thermodynamics and phase structure of an Abelian gauge 
theory at finite temperature that is rendered anisotropic by a finite density of charges which 
couple to a 2+1-dimensional Chern-Simons operator. The results obtained can be expected 
to remain qualitatively unchanged also at small coupling of its non- Abelian version. In 
maximally supersymmetric form the latter can be studied by gauge-gravity duality in the 
limit of infinite coupling and color number [32, 33]. 

Because of the resulting nontrivial phase structure, already the Abelian model studied 
here is certainly interesting in its own right. The physics motivation for considering an 
external density of charges with anisotropic coupling to gauge bosons is that it may provide 
a tractable toy model of anisotropics occurring in the (nonequilibrium) dynamics of quark- 
gluon plasma produced in heavy-ion collision, which is expected to be radically different 
at weak and strong coupling. 

In this paper, we have compared the resulting phase diagrams of this toy model in the 
weak and the strong coupling cases and found stark differences at high temperatures, but 
complete agreement in the limit of zero temperature (when expressed in terms of the scale 
do)- The differences are mainly due to the different signs of the free energy contributions 
of the Chern-Simons charge density a when a T. At weak coupling, a > introduces 
anisotropic mass terms for the gauge bosons which reduce the absolute value of the free 
energy compared to the isotropic Stefan-Boltzmann result. At strong coupling, nonzero a 
has the opposite effect. For the same reason, the pressure anisotropy at weak coupling is 
prolate in the high-temperature (low a) limit, while at strong coupling it is oblate. 

For T ;$ 0.3ao, the weak and the strong coupling limits lead to the same kind of 
instabilities of homogeneous phases in regions of their T-a phase diagrams corresponding 
to oblate pressure anisotropy. There are metastable and completely unstable regions where 
it is energetically favorable to have the charge density a redistributed along the axis of 
anisotropy so that a mixture of denser and empty (isotropic) regions arises in the form of 
stacks which are homogeneous in the transverse directions. The strong coupling results 
differ from the weak coupling ones in that the oblate phase extends to infinite temperature 
and so do the inhomogeneous phases appearing within that. At weak coupling, on the 
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other hand, the appearance of inhomogeneous phases is restricted to T and a both smaller 
than clq. Moreover, for moderate temperatures there are also mixed phases consisting of 
different prolate plasma; those phases have no counterpart in the strong coupling limit. 

Comparing to the situation of a weakly coupled plasma with momentum-space aniso- 
tropy we note that the instabilities of the present model are vaguely reminiscent of plasma 
instabilities. However, as we have seen, in the anisotropic Chern-Simons deformed theory 
they are not associated with unstable modes in the gauge fields. While this may point 
to limitations of this model as a model for the nonequilibrium situation of an anisotropic 
quark-gluon plasma, it also provides opportunities for further studies. As mentioned in the 
introduction, in the strong-coupling limit also several transport coefficients of relevance 
to heavy-ion phenomenology have been worked out already in the present model of an 
anisotropic plasma. It might be interesting to also obtain the corresponding results at weak 
coupling. Because of the absence of tachyonic modes in the gauge boson spectrum, such 
calculations are in principle feasible, while the hard-loop effective theory of a weakly coupled 
plasma with momentum-space anisotropy typically suffers from nonintegrable singularities 
[39, 42, 59]. 
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A. Computations 

In the presence of the source j(z) = a, the path integral in Eq. (2.1) can be carried out 
immediately with the result 

Z (T,a) = e-e™ ** K (A-HK)) ^ 
/det^ K (D^(K)) 



where the inverse photon and ghost propagators are given in momentum space by 

D~l{K) = (3 2 K 2 5 a p, (A.2) 

D7 l (K) = (3 2 [K 2 S tj ~ ae ljk k k ] , (A.3) 

D^{K) = D-}{K) = 0, (A.4) 

A^iK) = f3 2 K 2 . (A.5) 

The determinant of the inverse photon propagator over the Lorentz indices is given by 

\\D-l\\ = (ft 2 K 2 ) 1 ' 2,6 fi & K 2 [K 4 + (K 2 - kf)a 2] 



a 2 ± a 4 + 4a 2 k 2 
(f3 2 ) 4 - 2 '(K 2 ) 2 - 2 ' H[K 2 + * ll - I . (A.6) 



± 



It can immediately be seen that the ghost contribution cancels the contribution from the 
two a-independent photon modes. The logarithm of the partition function can thus be 
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written as 



where 



lnZ(T,a) = -^^^ln^fcg + ^^i + M^fc,,)))-^, (A.7) 



M^(fc,|) = A;} + 



a 2 ± Ja 4 + 4a 2 k 2 



or, equivalently, 



M ± (fc||) = Jk 2 + a 2 /4± 



(A.8) 



(A.9) 



Note that in Minkowski space (ko — > iuS) we have photon modes that for nonzero fey 
split into modes with uj\ > k 2 = k\ + A; 2 and ones with w 2 < k 2 . However, also the latter 
is non-tachyonic since w 2 > 0. 

Using 



a^ln^vrV+x 2 ) = 



1 



1 



4-7T 2 n 2 + x 2 2x 



1 



e x - 1 



= a x 2 [x + 21n(l-e- ;E )] 

oo 

^ ln(4vr 2 n 2 + x 2 ) = x + 2 In (l - e~ x ) + constant, 



(A.10) 
(AH) 



we get for the Helmholtz free energy density / = — T/VlnZ (with the constant absorbed 
into the cosmological constant) 



f(T, a) = O + T W \Uu± + In (l - e"^±) 

_l_ */ k L 



where 



w| = fcl + M|(fc||). 

A.l T = contribution 

The T = contribution to the free energy density is given by 
1 



(A.12) 



(A.13) 



/(0,o) = fi+- / [W+ (fc_L , fc[| ) + 6LF_ (fej. , fc|| )] 



(A.14) 



4ir 



fen JO 



dk±k± 



kl + M 2 + m) + Jkl + M 2 _{k\\) 



where we have introduced a UV regulator for the integration over the transverse momenta 
and anticipated that the zero-loop contribution (which is just fi) will need to be renormal- 
ized. The integral over the transverse momentum can be carried out easily, giving 

r°° " M±(k\\) 



± dk ± k ± ^k 2 ± + Ml{k\\) = dkk yj 'k 2 + Ml{k\\) - k - 



2k 



+ dk [k 2 + -M 2 ± {k\\) +C(1/A 



(A.15) 



: M(k\\) + ^Ai + -M 2 ± (fc||)A ± + 0(1/Aj 
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Inserting this into /(0, a) above, dropping terms that vanish as the regulators are removed, 
we get 

/(0 ' a) = "l27rl [ M +W + M -W]+^ I (2fcJ + a 2 )+^(A ± ,e) 

(A.16) 



12 _ , [Ml(k\\) + Mi(k\\)]+n B (Aj_,e). 



Extracting the large-fen behaviour, the first term can be written as 



£ L M ^ = I f 

± Jk n 71 Jo 



dkw 



M%{k\\) + Af£(fc||) - 2fef - -a 2 fej| 

+ 0(e) 



15a 4 



64jfc 2 + a 2 



,2 9 2 , 15a 4 
" 4 " 64jfe 2 + a 2 



+ 



vr 128tt V e 



ca 4 15a 4 ( \ , a 2 . 

• . . f 1" A , ) K?(r). 



(A.17) 



where 



jHda; (x 2 + ^(1 + Vl + 4x 2 )^) 7 + [x 2 + ^(1 - y/l + 4x 2 ) 



\ 3/2 



-2x d - -x - 



15 



4 64V1 + a; 2 
The free energy density at T = is then 

ca 4 5a 4 fl 

/(0,a) = 



0.29136... 



(A.18) 



12** 5i2*n* ln Fj + ^' e) 

4 c 4 

co* 5a^ , a 



12vr 2 256tt 2 A 
where f2(A) is the renormalized cosmological constant in the MS-scheme. 

A. 2 Finite T contribution 

The finite-T contribution is given by 

f(T,a)-f(0,a) = Tj2 [ In (l - e"**) . 



(A.19) 



(A.20) 



The integral over the transverse momenta can be carried out immediately using standard 
integrals. Noting that 



2tJ^ In ( 1 - e~^ ± 
dk 2 t ' 



1 1 



oj± e 



fa± - 1 ' 



(A.21) 
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we find 

d 2 k { , / \ 3T 3 



A* i / \ 3T d 

_^ln(l-e-^) = — MM ± (fc||)/T,0), (A.22) 

where the function h n (y,r) is defined by 

1 f°° x n ~ x 1 

h n (y,r) = —, . / c?x — . — . . (A. 23) 

nVy ' 7 r(n)7 y/x 2 + y 2 e V^-ry _ i 

Their properties are discussed in detail in Ref. [60]. In particular, 

—h n+1 (y,0) = -^-h n - 1 (y,0), (A.24) 
ay n 

h 2 (y,0) = -ln(l-e-»), (A.25) 
Mtf, 0) = ^ [y Li 2 (e-v) + Li 3 (e-w)] . (A.26) 
We now get for the integral over the longitudinal momentum, 

-^h 4 (M±/T, 0) = - / dku h 4 (M±/T, 0) 

= ^ dk^k\\M ± d ^h 2 (M ± /T,Q) (A.27) 



- rdxx 2 [l± - V llnfl-e-V^^ 



3tt 7o V \J x 2 + 

where we have denoted y = a/ (2T). The finite-T contribution is thus given by the integral 

/<r,„)-/(0,.) = ^jT^^i^^J^l-.-^*) . (A.28) 

A. 2.1 High- and low-T expansions 

Even though we shall study the thermodynamic properties of this system numerically, it 
is instructive to calculate the analytic high- and low-T approximations to the free energy 
as well. As derived above, the finite-T contribution is given by 

/(T, a) - /(0, a) = --2 - [h s (y, 1) + h(y, -1)] + ^I(y, 1), (A.29) 

where 

7( y , r ) = dx ^ V f = [in (l - e -V^+i^-n/^ _ m (i _ e -V* 2 +2/ 2 +n/jj . (A.30) 
We immediately find that 

l(y,0)=0, (A.31) 
(9/ 

— = 2y 2 [h 3 (y,r) + h 3 (y, -r)] (A.32) 
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and thus 

I(y, 1) = 2y 2 [ dr [h 3 (y, r) + h 3 (y, -r)] . (A.33) 
J o 

The high-T (small-?/) expansions of h n (y,r) + h n (y, —r) have been computed in Ref. [60], 
with the results 

h(y, 1) + h 5 (y, -1) = — + — + — ( In + 7jB + — ^ + 0(y 6 ), (A.34) 
M^ r) + h 3 (y, -r) = y - ^ - ^ (in J- + 7 £ - ^ + r 2 ) + 0(y 4 ) (A.35) 
Putting everything together, and setting y = a/(2T) we get 

/(T,„)-/(0, o) = -_+_-_-_(,*— +7E --) +0 („<T-) (A.36) 
for the high-T expansion. On the other hand, at low temperatures (high-y limit) we have 
Ml/. 1) + My, -1) = (8C(5/2)y + 15C(7/2) + 0(1/ yj) (A.37) 

Ml/.O + Ml/,-') = ^/^fLi3/ 2 (e (r - 1)s/ ) + Li 3/2 (e-( r+1 ^) 



2V 2 
3 

/* 
8V"2 



+ A ( L i 5/2 (e^) + Li 5/2 (e-^)) +0(l/y 2 ) 



(A.38) 

^ 1) = (8C(5/2)y + 3C(7/2) + 0(l/y)) (A.39) 



Combining the results gives us 



f(T, a) - /(0, a) = -^?« 1/2 T 7 / 2 + O(a^) (A.40) 



. 3C(7/2) „l/2 T 7/2 + ro ^-l/2 

in the low-T limit. 
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